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We present a filter for detecting gravitational wave signals from burst sources. This filter requires 
only minimal advance knowledge of the expected signal: i.e. the signal's frequency band and time 
duration. It consists of a threshold on the total power in the data stream in the specified signal band 
during the specified time. This filter is optimal (in the Neyman-Pearson sense) for signal searches 
where only this minimal information is available. 



I. INTRODUCTION 

Currently, the best understood and most highly devel- 
oped technique for detecting gravitational waves with in- 
terferometric detectors is matched filtering. This tech- 
nique is optimal if the waveform to be detected is known 
in advance. There are, however, potentially impor- 
tant sources of gravitational radiation that are not well 
enough modeled to obtain reliable waveforms. Included 
in this category are binary black hole mergers, which have 
been discussed in some detail by Flanagan and Hughes 
(FH), and supernovae 

In order to detect poorly modeled sources, new tech- 
niques must be developed. Clearly, these techniques must 
perform well with incomplete prior knowledge of the ex- 
pected signal. A number of techniques are under active 
investigation ||^. 

The purpose of this article is to discuss one such tech- 
nique, the power filter . This filter only requires prior 
knowledge of the duration and frequency band of the sig- 
nal. It is therefore well suited to detecting black hole 
merger signals, for which FH have estimated these pa- 
rameters. Furthermore, we shall show in the following 
sections that this filter is optimal in the sense that it 
gives the highest probability of correctly detecting a sig- 
nal for a given false alarm probability. Our treatment 
here is cursory; a more comprehensive description of the 
filter is in preparation H. 



II. THE POWER STATISTIC 

Consider the output h{t) of the interferometric gravi- 
tational wave detector. It is sampled at a finite rate 
1/At to produce a time series hj = h{jAt), where j is 
an integer. A segment of TV samples defines a vector 



h={h 



■J' ■ 



h 



This vector can be written as 



(2.1) 



where n is detector noise and s is a (possibly absent) 
signal. 



The noise is assumed to be stochastic. The vectors 
n and h are therefore described statistically. Statistical 
fiuctuations lead to two types of errors in detecting a 
signal: false alarms, in which signals are detected when 
none are present, and false dismissals, in which signals are 
not detected when present. An optimal filter is defined 
to be one which minimizes false dismissals for a given 
false alarm rate. 

Neyman and Pearson have shown that an optimal filter 
is one for which the likelihood ratio A is compared to a 
threshold [plJq| . The likelihood ratio is defined to be 



A[h] ^ / P[s] 



p[h|s] 
p[h|0]' 



(2.2) 



where p[h|s] (p[h|0]) is the probability of obtaining h 
given that a signal s is present (absent) and ^[s] is a 
measure on the space of signals. 

The quantities p[h|s] and p[h|0] depend on the statis- 
tical properties of the noise. For convenience, we assume 
in this article that the noise is stationary and Gaussian^. 
We can therefore write the probability distribution for 
the noise as 



p[n] — C exp 



n • n 



(2.3) 



where C is a constant prefactor and n ■ n an inner prod- 
uct, both of which are determined by the autocorrellation 
matrix of the noise. When a signal is present (abs ent) we 
have n = h — s (n = h) , and can easily use Eq. ( |2.3| ) to 
find the integrand in Eq. (2.2) 



p[h I 0] 



= exp 



(s.h)-i(s. 



(2.4) 



The measure I?[s] in Eq. (2.2) reflects our prior knowl- 
edge about the signal. For those signal parameters about 
which we have no prior knowledge, we choose a measure 
which reflects our ignorance. 



However, other types of noise can also be considered 
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Consider now the case where we know the time win- 
dow and frequency band in which a signal occurs. The 
measure V[s] restricts the integral in Eq. ( ^.2[ ) to the 
projection h|| of h into the space of vectors with the give 
window and frequency band. Introducing the notation 
= s • s, £ = h|| • h|| and s • hy = AS^^^ cos 9, we 
rewrite Eq. (^.2[) as 



A[h] 



V[9,A] exp 



A£'^^ cos9--A^ 
2 



(2.5) 



Since we claim no prior knowledge of 9, a suitable mea- 
sure over 9 is uniform over all possible angles between h|| 
and s (i.e. over an J\f sphere, where J\f is the dimension 
of the space of vectors with the required duration and 
frequency band), which reflects our lack of knowledge. 

While a suitable measure over the signal amplitude A 
can also be deduced, it is unnecessary here. Instead, one 
constructs a locally optimal statistic [^, Aioc[h], which 
is appropriate in the limit of weak signals. To construct 
this statistic, one expands the likelihood ratio (2.5) in a 
Taylor series about A — 0. The statistic is simply the first 
non- vanishing coefhcient (excluding the A° coefficient) in 
the expansion. Expanding (|]^) and integrating over 9 
we get 



Aioc[h] oa £ + terms independent of h. 



(2.6) 



The terms independent of h clearly do not discriminate 
between the presence and absence of a signal. Thus, the 
optimal statistic for detecting a signal of known dura- 
tion and frequency band is simply the total power in the 
detector output over that time and band. 



III. OPERATING CHARACTERISTICS 

In the previous section we determined the optimal statis- 
tic £ for signals of known duration and frequency band in 
stationary Gaussian noise. We construct the optimal fil- 
ter from this statistic via a threshold decision rule. That 
is, we calculate at what value £* of the statistic we incur 
the largest acceptable false alarm probability. We then 
compare values of £ calculated from our data with £*. A 
signal is said to have been detected if £ > £* . 

For Gaussian noise, false alarm and false dismissal 
probabilities can be calculated analytically up to quadra- 
ture. If no signal is present, £ is just the sum of the 
squares oiV = 2x6txSf independent Gaussian random 
variables. Thus £ has a distribution with V degrees 
of freedom, and the false alarm probability for a value £* 
is just 



P{£>£* \A^ 0) 



r(y/2^fV2) 



(3.1) 



where T{a, x) = J°° e *'t°' ^dt is the incomplete Gamma 
function. 



If a signal of amplitude A is present, £ is distributed 
as a weighted sum of probability distributions, 

p{£ \V^A) = Y, e-^V^(AV2)" e-^/2(£:/2)"+W2-i 



n=0 



r(n + y/2) 



(3.2) 



This is the non-central probability distribution. The 
false dismissal probability is given by 

P{£ <£* \A)^ [ p{£\ V, A) d£ (3.3) 
This probability can be integrated numerically. 

IV. SUMMARY AND DISCUSSION 

We have presented here a power filter to search for gravi- 
tational wave signals from burst sources in interferomet- 
ric data. The filter is designed to look for signals of 
known duration and frequency bandwidth; when this is 
the only available information, and the noise is stationary 
and Gaussian, the power filter is optimal. Moreover, this 
filter is locally optimal for a wide class of non-Gaussian 
noise, thus making it a useful tool to analyze real inter- 
ferometer data. 

One shortcoming of the power filter is its inability to 
distinguish between gravitational wave signals, and spu- 
rious instrumental artifacts which produce time and band 
limited signals. This is mitigated by using multiple in- 
struments to detect gravitational wave bursts. As an 
added benefit, bursts identified as noise can then be used 
for detector characterization. The extension of the power 
filter to multiple instruments will appear in an upcoming 
article Q ; this article also contains a more complete dis- 
cussion of the filter, including implementation strategies 
and a comparison to matched filtering. 

In conclusion, we think that the power filter provides a 
useful tool for gravitational wave data analysis. It should 
play a significant role in detector characterization for sin- 
gle interferometers, and should form the basic building 
block in an hierarchical detection strategy using multiple 
interferometers. 
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